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We use the Markov Chain Monte Carlo method to investigate a global constraints on the modified 
Chaplygin gas (MCG) model as the unification of dark matter and dark energy from the latest 
observational data: the Union2 dataset of type supernovae la (SNIa), the observational Hubble data 
(OHD), the cluster X-ray gas mass fraction, the baryon acoustic oscillation (BAO), and the cosmic 
microwave background (CMB) data. In a fiat universe, the constraint results for MCG model are, 
n,h^ = 0.02263l°:»i»^ (la) (2a), = 0.7788l°:g^,l (2-), a = 0.1079i«il^ 

(1-) t'.fJA (2a), B = 0.00189i°:S (1-) t^'.^s?^ (2-), and Ho = 70.71ltt.}^^ (1-) (2-)- 



PACS numbers: 98.80.-k 
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I. Introduction 



^ . Recently, mounting cosmic observations suggest that the expansion of present universe is speeding up rather than 
I slowing down 1] . And it indicates that baryon matter component is about 5% of the total energy density, and about 
■ 95% of the energy density in the universe is invisible, including dark matter (DM) and dark energy (DE). In addition, 
it is shown that DE takes up about two-thirds of the total energy density from cosmic observations. In theory many 
' kinds of DE models [2] have already been constructed to explore the DE properties. 

Chaplygin gas (CG) and its generalized model have been widely studied for interpreting the accelerating universe 
01 . The CG model can be obtained from the string Nambu-Goto action in the light cone coordinate For 
. generalized Chaplygin gas (GCG), it emerges as a effective fluid of a generalized d-brane in a (d -I- 1, 1) space time, 

H ■ . . n . . 

d • and its action can be written as a generalized Born-Infeld form 6] . Considering that the application of string theory in 
principle is in very high energy when the quantum effects is important in early universe [Si] , the quantum cosmological 
studis of the CG and the GCG has been well investigated in Ref. Q and Q- In addition, one knows that the most 
attractive property for these models is, two unknown dark sections in universe-dark energy and dark matter can be 
unified by using an exotic equation of state. It is worthwhile to study the unified models of dark sections for other 
generalization of CG. 
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A simple and popular generalization relative to the GCG model is that it is extended to a form by adding a 
barotropic term, referred to as the modified Chaplyein gas (MCG) 3- The correspondences between the MCG and 
the ordinary scalar field [8] , the tachyon theory "3, Q , and the holographic dark energy density {lo| of the universe 
have been studied. In principle, MGG represents the evolution of the universe starting from the radiation era to the 



era dominated by the cosmological constant lll-ll3l|. In addition, MCG is also applied to inflation theory [U]. For 



the more discussion on MCG model, please see Refs. 
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ITj . In this paper, we use the Markov Chain Monte Carlo 



(MCMG) technique to constrain this more general unified candidate, MCG. The used observational data include, the 
Union2 data of type la supernovae (SNIa) 18 1, the observational Hubble data (OHD) jl^, the cluster X-ray gas mass 



fraction [20|, the measurement results of baryon acoustic oscillation (BAG) from Sloan Digital Sky Survey (SDSS) 



and Two Degree Field Galaxy Redshift Survey (2dFGRS) [21 1 



(CMB) data from seven-year WMAP [23 1 



221, the current cosmic microwave background 



II. Modified Chaplygin gas model 

We briefiy introduce the GCG model as the unification of dark matter and dark energy at first. The energy density 
p and pressure p in this model are related by the equation of state (EOS) 

PGCG = a ' (1) 

Pgcg 

where A and a are parameters in the model. When a = 1, it is reduced to the CG scenario. Using Eq. ([IJ one has a 
solution of energy density for the GCG fluid 

1 — A 1 

PGCG = POGCG [As + ^^i+a) ] ~ ' (2) 

where As = i+a, ■ From Eq. ([2]) one can see that the GCG fluid behaves as a dust-like matter at early time and as 

PoGCG 

a dynamical cosmological constant at late epoch [24], then the GCG can be interpreted as an entangled mixture of 
dark matter and dark energy. The dual role of the GCG fluid is at the heart of the interesting property of this model. 

For MCG model, from phenomenological view point it is interesting and can be motivated by the brane world 
interpretation [25]. It is characterized by a more general EOS, 

Pmcg = BpMcc 5 , (3) 

Pmcg 

which looks like that of two fluids, one obeying a perfect EOS p — Bp and the other being the GCG [sj. Where A, B, 
and a are parameters in the model, pmcg f^nd pmcg £^re energy density and pressure of the MCG fluid. From Eq. 
([3]), it is easy to see that the EOS of MGG reduces to the GCG scenario if i? = 0, and reduces to the perfect fluid if 
A = 0. In addition, with B = Q and a = it reduces to AGDM model. In the space-time geometry described by the 
non-flat Friedmann-Robcrtson- Walker (FRW) metric 

ds^ = -dt^ + a^{t)[- — + r^ide^ + sin^ 0d^% (4) 

the EOS ((21) leads, after inserted into the relativistic energy conservation equation, to an energy density of MCG 
evolving as 

Pmcg = Pomcg [Bs + ^3(i+B)(i+a) ] 
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for B 7^ -1, where Bs — - — "^i i „ — , a is the scale factor of universe which is related to the redshift by, a = -r^. 

(l + -«)P0AfCG 

Considering that MCG fluid plays the role of a mixture of dark matter and dark energy, and assuming that the 
universe is filled with three components: the MCG, the baryon matter, and the radiation component, one can express 
the dimensionless Hubble parameter E as 

E=^ = Jil-^i-^r- nk)[Bs + (1 - S,)a-3(i+B)(i+")]TT^ + r!ba-3 + r!,a-4 + 17fea-2, (6) 

^0 

where H is the Hubble parameter, with its current value Hq = 100ft- km s~^Mpc~-'^, fif,, fir, and ftk denote the 
dimensionless baryon matter, radiation, and curvature density, respectively. 



III. Constraint result and conclusion 

Next we apply the current observed data to constrain the MCG model. The constraint method and the used data: 
Union2 SNIa, OHD, CBF, BAO, and CMB data are presented in Appendix A. In our calculations the total likelihood 
function is written as L oc e^^ Z^, with the total equaling 

X = XSNIa + XOHD + XCBF + XbAO + XcMBi \^) 

where the separate likelihoods of SNIa, OHD, CBF, BAO and CMB are given by Eqs (jMI, (EHl), (EUl, (|X24t and 
(|M8l) . From the expressions of XcB-F' Xbao ^^'^ Xcmb^ one can see that they are related with the matter density 
flm- For MCG model, since it is considered as the unification of dark matter and dark energy, we do not have dark 
matter in this model. So, the matter density is not explicitly included in the background equation (jG]). According the 
Eq. ([5]), by considering the universe is dominated by the matter component at early time (a <C 1), i.e., relative to the 
dark matter density the dark energy density is neglectable, one can get an effective expression of the current matter 
density, i7„j = ilb + (1 — il^ — ilr — i^fe)(l — ^s)^^ • This expression of i7„i, is an estimate of the "matter" component 
of the MCG fluid with the baryon density. Thus, in the CBF, BAO and CMB constraints, we take this expression of 

Next, we perform a global fitting on determining the MCG model parameters using the Markov Chain Monte Carlo 
(MCMC) approach. In our joint analysis, the MCMC code is based on the publicly available CosmoMC package [2^ 



and the modified CosmoMC package [20|, |27|, |28| . The latter package is about the constraint code of X-ray cluster 
gas mass fraction, with including additional 7 free parameters {K,ri,^,bQ, a^, sq, as). In these packages, they have 
been modified to include the new model parameters B^, a and B. In the calculation the baryon matter density is taken 
to be varied with a tophat prior: ilbh^ E [0.005, 0.1]. In addition, for the MCMC calculation on the MCG model, we 
run 8 independent chains, and to get the converged results we test the convergence of the chains by typically getting 
i? — 1 to be less than 0.03. 

In FIG. [U we show a one-dimensional probability distribution of each parameter and two-dimensional plots for 
parameters between each other in the flat MCG model. According to the figure, the constraint results on the best 
fit values of parameters with la and 2a confidence levels are listed in Table H] In Ref. ^13^ it is shown that for the 
constraint on MCG model parameters obtained from the sound speed, < < c^, the parameters are restricted to 
< (B -1-1) (a -1-1) < 2. According to this constraint on MCG model parameters one can see that the MCG fiuid has a 
well sound speed. Considering Ref. [l^ 1 where the model parameter B is constrained by using the location of the peak 
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Parameters Best fit values 





U.UZZOO_QQoig2_o. 00195 


Bs 


n 77CC+0. 0736+0. 0918 
U.( 'oO_o. 0723-0. 0904 


a 


n 1070+0.3397+0.4678 
U.lU(y_g 2539-0.2911 


B 


n nni SQ+o oosss+o. 00660 

U.UU10y_Q Qg7gg_g QQgjg 


Ho 


7n 71 1 +4.188+5.281 
/U. 1 11_3,142-4.149 



xLuixlir^/dof) 596.147 (0.9709) 



TABLE I: The data fitting results of the MCG model parameters with la and 2a confidence levels. 



of the CMB spectrum, —0.35 ^ B < 0.025, it is easy to see that our result is more stringent. In addition, considering 
Ref. [29I , where the values of model parameters are analyzed against the matter power spectrum observational data, 
it is obtained that a very stringent constraint exists on B, which is consistent with our result, i.e., it seems that the 
MCG model is viable only for very special cases, and it tends to reduce to the GCG scenario. 




D.O22 0,024 






FIG. 1: The 2-D contours with la, 2a confidence levels and 1-D marginalized distribution of Qbh^, -Bs, a, B, and Ho in the 
flat MCG model. 

In addition, we consider a more general non-flat background geometry. For using above combined observational 
data to constrain the MCG model, since the space curvature fi^ is near to zero, the constraint results of the model 
parameters obtained from the non-flat universe are similar to the cases of the flat universe. For simplicity, we do not 



list the constraint results for this case. With replacing the Union2 SNIa data-'^ with the 397 Constitution data^ [30| in 

the above combined constraint, we obtain the constraint results of MCG model parameters, flk — — 0.000844tooi5i9i 
+0.014536 /r,„\ n n 7t;/n +0.0941 +0.1092 fn\ ^ n nnfio+0.3837 +0.4401 fr,\ E) n nm Qs+0 00817 

(1^) 10.017862 (2c^)' = 0.7541_o 0892 (lo-) _o.0965 (2cr), a = 0.0082_o 3318 (1^) -0.2433 (^cr), B = 0.00138_o.O0738 

(Icr) to:oo869 (2^^), with xl^^n = 520.118 {xlmn/dof = 1.148). It sccms that for this case, the MCG model tends to 
reduce to the flat cosmic concordance model, ACDM. 

Acknowledgments The research work is supported by the National Natural Science Foundation (Grant No. 
10875056), NSF (10703001) and NSF (No.ll005088) of P.R. China. 



Appendix A: Current observational data and cosmological constraints 

In this part we introduce the cosmological constraint methods and the current observed data used in this paper. 

1. Type la supernovae 

SNIa behave as the excellent standard candles, so they can be used to directly measure the expansion rate of the 
universe from the high redshift to the present time. For using SNIa data, theoretical dark-energy model parameters 
are determined by minimizing the quantity [sil 



XsNIa (MO , ) = 2^ 2 r~\ ' '^^^> 

where N — 557 for Union2 dataset, which is the largest SNIa sample by far; ps denotes the model parameters; cr^^j^^ (zi) 
are errors; fJ.obsizi) — mobs{zi) — M, is the observed value of distance modulus of SNIa at Zi and can be given by the 
SNIa dataset; ^th is the theoretical distance modulus, which is related to the apparent magnitude of SNIa at peak 
brightness m and the absolute magnitude M, 

^ithiz) = mth{z) - M = 5logio{DLiz)) + ^o- (A2) 

Here, the Hubble free luminosity distance 

DLiz) = HodUz) = Sil±^,,nn[^\ f f",^^ ], (A3) 
V l^'fcl "'0 \z ,Ps) 

and 

TT-l 

Mo - 5Zo5io(-^) + 25 = 42.38 - blogioh, (A4) 
Mpc 

with h being a re-normalized quantity, which is given by Hq = lOOh km s^^Mpc^^. It should be noted that /io is 
independent of the data and the dataset. By expanding the of Eq. (jAip relative to the nuisance parameter hq, 
the minimization with respect to /xg can be made trivially |32| 

xlNiaiPs) = A{ps) - 2ixaB{ps) + flic, (A5) 



^ The Union2 SNIa data are obtained, by adding new datapoints (including the high redshift SNIa) to the Union SNIa data, making a 

number of refinements to the Union analysis chain, refitting all light curves with the SALT2 fitter. 
^ The 397 Constitution data are obtained by adding 90 SNIa from CfA3 sample to 307 SNIa Union sample. CfA3 sample are all from the 

low-redshift SNIa, z < 0.08, and these 90 SNIa are calculated with using the same Union cuts. 
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where 



N 



1=1 



(A6) 



TV 



^J■obs{Zi) - 1-J.thiZi] fJ-O ^ 0,Ps) 



(A7) 



(A8) 



Obviously, according to Eq. (jASP XsNia has a minimum for — B/C. Thus, the expression of x for SNIa constraint 



can be written as 



xInM = MPs) - Bips)yc. 



(A9) 



Since x'sNia min — XsNia min ^nd XsNia independent of the nuisance parameter /io, one usually utilize the expres- 
sion (IA9|) to displace (jAip to perform the likelihood analysis for the SNIa constraint. For minimizing XsNiaiPs^ B/C) 
to constrain cosmological model, it is equivalent to maximizing the likelihood 

-X^iPs)^ 



L{ps) oc exp[ 



(AlO) 



Observational Hubble data 



The observational Hubble data [33| are based on differential ages of the galaxies. In 3J], Jimenez et al. obtained 



35l |. and used it to constrain the 



an independent estimate for the Hubble parameter using the method developed in 
cosmological models. The Hubble parameter depending on the differential ages as a function of redshift z can be 
written in the form of 



H{z) = - 



1 dz 



(All) 



1 + z dt 

So, once dz/dt is known, H{z) is obtained directly. By using the differential ages of passively-evolving galaxies from 



the Gemini Deep Deep Survey (GDDS) [36] and archival data |37H42| | . Simon et al. obtained several values of H{z) 



at different redshift 



19[. The twelve observational Hubble data (redshift interval < z < 1.8) from Refs. 43|, |4J| are 



listed in Table |TI1 In addition, in 



45| the authors take the BAO scale as a standard ruler in the radial direction, and 



z 





0.1 


0.17 


0.27 


0.4 


0.48 


0.8J 


5 0.9 


1.30 


1.43 


1.53 


1.75 


H{z) (km s-^Mpc"^) 


74.2 


69 


83 


77 


95 


97 


90 


117 


168 


177 


140 


202 


1(T uncertainty 


±3.6 


±12 


±8 


±14 


±17 


±60 


±40 ±23 


±17 


±18 


±14 


±40 



TABLE II: The observational H{z) data [4a,|44 
obtain three additional data: (z = 0.24) = 79.69 ± 2.32, H{z^ 0.34) = 83.8 ± 2.96, and i?(z = 0.43) = 86.45 ±3.27. 
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The best fit values of tfie model parameters from observational Hubble data are determined by minimizing 46 1 



where Htu is the predicted value of the Hubble parameter, Hobs is the observed value, a{zi) is the standard deviation 
measurement uncertainty, and the summation is over the 15 observational Hubble data points at redshifts Zi. 



3. The X-ray gas mass fraction 

The X-ray gas mass fraction, /gas, is defined as the ratio of the X-ray gas mass to the total mass of a cluster, which 
is approximately independent on the redshift for the hot {kT > ^keV), dynamically relaxed clusters at the radii larger 
than the innermost core r25oo- As inspected in [2^, the ACDM model is very favored and has been chosen as the 
reference cosmology. The model fitted to the reference ACDM data is presented as [2^ 



fACDMf^ _ KA-fb{z) ( nt 

J nas \^) 



Da{z) 



(A13) 



where D')^'^^'^ [z) and Da{z) denote respectively the proper angular diameter distance in the ACDM reference cos- 
mology and the current constraint model. A is the angular correction factor, which is caused by the change in angle 
for the current test model ^2500 in comparison with that of the reference cosmology ^2oto'*^ ■ 



\ 02500 J ~ \ [H{z)Da{z)]^^^^ 



A^ l^ ^ , (A14) 



here, the index r/ is the slope of the f gas{i' / ^2^00) data within the radius r2500 7 with the best-fit average value 
77 — 0.214 ± 0.022 ^]. And the proper (not comoving) angular diameter distance is given by 

Da{z) = - ^^sinn[vW/ tTT^]- (A15) 

{l + z)^/\Sh\ Jo H{z';ps) 

which is related with dL{z) by 

d-Liz) 



Da{z) 



(1-f z)2' 

where sinn{yj\nu\x), respectively, denotes s\iv{^J \Vlk\x) , ^J\^k\x, sm)[Y{^J\D,k\x) for Vtk < 0, fife = and Qk > 0. 

In equation (|A13P , the parameter 7 denotes permissible departures from the assumption of hydrostatic equilibrium, 
due to non-thermal pressure support; the bias factor b{z) = bo{l + ai,z) accounts for uncertainties in the cluster 
depletion factor; s{z) = so(l -I- Ugz) accounts for uncertainties of the baryonic mass fraction in stars and a Gaussian 
prior for sq is employed, with sq = (0.16 ± 0.05)/iy5^ 2^; the factor K is used to describe the combined effects of the 
residual uncertainties, such as the instrumental calibration and certain X-ray modelling issues, and a Gaussian prior 
for the 'calibration' factor is considered by K ^ 1.0 ± 0.1 [2(|. 
Following the method in Ref. 



20 



471 and adopting the updated 42 observational fgas data in Ref. '20], the best 



fit values of the model parameters for the X-ray gas mass fraction analysis are determined by minimizing, 

2 _^ [f^a's'''H^^)-f9as{z^)? , (^0-0.16)^ , (i^^l-O)^ , (r? ~ 0.214)^ . ^ ^ ^, 

^CBF-l^ ^^j^^ + 0.00162 + 0.012 + 0.0222 ' iA16j 

i Jgas 



where o-fg^^{zi) is the statistical uncertainties (Table 3 of [2(|). As pointed out in t20|, the acquiescent systematic 
uncertainties have been considered according to the parameters i.e. 77, b{z), s{z) and K. 
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4. Baryon acoustic oscillation 

The baryon acoustic oscillations are detected in the clustering of the combined 2dFGRS and SDSS main galaxy 
samples, which measure the distance-redshift relation at zbao = 0-2 and zsao — 0.35. The observed scale of the 
BAO calculated from these samples, are analyzed using estimates of the correlated errors to constrain the form of the 



distance measure Dv{z) 2lA, 3 



Dv{z) = 



{l + zfD\{z) 



H{z) 



1/3 



(A17) 



The peak positions of the BAO depend on the ratio of Dv{z) to the sound horizon size at the drag epoch (where 
baryons were released from photons) z^, which can be obtained by using a fitting formula [4^: 

1291(l]™/^2)-0.419 



1 + 0.659(O„/i2) 



0.828 



2^^)2l 



with 



hi = 0.313(r!™/i2)-o.4i9j^ ^ 0.607(17™/i2)0.674^ 
&2 = 0.238(17™/i2)0-223_ 



(A18) 

(A19) 
(A20) 



In this paper, we use the data of rs{zd) / Dv{z) extracted from the Sloan Digitial Sky Survey (SDSS) and the Two 
Degree Field Galaxy Redshift Survey (2dFGRS) [48|, which are listed in Table Hill where rs{z) is the comoving sound 
horizon size 



, Csdt f Csda f°° Cs 
r, z =c / — c / — — c dz — -— 



1/(1 + 2) 



da 



where Cs is the sound speed of the photon— baryon fluid 



V3J0 a?H{a)^l + (3r2f,/(4r2.^)a) ' 



50 



c- = 3+lxM£l.3 + lx(g^)a, 

3 P7(z) 6 



(A21) 



(A22) 



and here = 2.469 x IQ-^h'^ for Tcmb = 2.72hK 



z 




0.2 


0.1905 ±0.0061 


0.35 


0.1097 ±0.0036 



TABLE III: The observational rs(zd)/7?v(z) data 



Using the data of BAO in Table Hill and the inverse covariance matrix V ^ in 22 1 

30124.1 -17226.9 



-YlTl%.'i 86976.6 



(A23) 



tJ^e XsAobs) i*^ gi'^en as 



(A24) 
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7- year maximum likelihood 


error, a 


Ia{z,) 


302.09 


0.76 


R{z,) 


1.725 


0.018 


z* 


1091.3 


0.91 



TABLE IV: The values of Ia{z*), -R(z*), and from 7-year WMAP data. 



where X is a column vector formed from the values of theory minus the corresponding observational data, with 

, n°^,n''l^ - 0-1905 , 
X = , (A25) 



Dv(0.35} 



-0.1097 



and X* denotes its transpose. 



5. Cosmic microwave baclcground 



The CMB shift parameter R is provided by 







R{z,) = ^n„,H^{l + z,)Da{z,)/c, (A26) 



here, the redshift (the decoupling epoch of photons) is obtained using the fitting function 53 1 

= 1048 [i + omi24{ni,h^)-° '^^^] [1 + giin^hy^] , 

where the functions gi and g2 read 

51 = 0.0783(r!b/l2)-0.238 (1^ 39 5(^^/^2)0.763)-!^ 



52 = 0.560 (1 + 21. 1(r2b/i 



2\1.81 



-1 



In addition, the acoustic scale is related to a distance ratio, DA{z)/rs{z), and at decoupling epoch it is defined as 

U^il + z,.)^^^j^, (A27) 

where Eq. (|A27|) arises a factor 1 + z*, because DaIz) is the proper (physical) angular diameter distance, whereas rs(z*) 
is the comoving sound horizon. Using the data of IatRtZ^, in [23] and their covariance matrix of [^^^(z*), _R(z*), z*] 
(please see table HVl and |V| . we can calculate the likelihood L as Xcmb ~ ~21nL: 

xImb - Ad,[Cov-\d,,dj)[Ad,]% (A28) 

where Adi = di — df""^"" is a row vector, and di = {Iai R,z^,). 
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